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Phenomenological models of the dipole cross section that enters in the description of for instance 
deep inelastic scattering at very high energies have had considerable success in describing the avail- 
able small- a; data in both the saturation region and the so-called extended geometric scaling (EGS) 
region. We investigate to what extent such models are compatible with the numerical solutions of 
the Balitsky-Kovchegov (BK) equation which is expected to describe the nonlinear evolution in x 
of the dipole cross section in these momentum regions. We find that in the EGS region the BK 
equation yields results that are qualitatively different from those of phenomenological studies. In 
particular, geometric scaling around the saturation scale is only obtained at asymptotic rapidities. 
We find that in this limit, the anomalous dimension y(r, x) of phenomenological models approaches 
a limiting function that is universal for a large range of initial conditions. At the saturation scale, 
this function equals approximately 0.44, in contrast to the value 0.628 commonly used in the models. 
We further investigate the dependence of these results on the starting distribution, the small-?- limit 
of the anomalous dimension for fixed rapidities and the ^-dependence of the saturation scale. 
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> 1 I. INTRODUCTION 

G\ '. 

At very high energy, deep inelastic scattering (DIS) of a virtual photon off a proton can be described as the 
scattering of a color dipole off small- a; partons, predominantly gluons, in the proton The evolution in x of 
the corresponding dipole cross section has been the subject of many investigations. The earliest and best-known 
£--«. . equation describing x-evolution at high energy, i.e. small values of x, is the BFKL equation 0,0]. It is a linear 
evolution equation based on resummation of logarithms in 1/x, which arc dominant at small x. The BFKL 
equation predicts an exponential growth of the dipole cross section as x decreases, which potentially violates 
unitarity. Attempts to solve this problem have resulted in the formulation of a variety of nonlinear evolution 
equations, the earliest example of which is the GLR equation [J, [|| for the gluon distribution. More recently, a 
nonlinear evolution equation for the dipole cross section has been obtained, the Balitsky-Kovchegov (BK) equation 
[1,0]. As a consequence of the nonlinearity, the dipole cross section saturates with decreasing x, thereby offering 
a resolution to the unitarity problem arising from the BFKL equation. 
. £h ' The momentum scale at which nonlinearity becomes important and saturation sets in is called the saturation 

scale Q s . This scale grows as x decreases, which means that if Q s becomes larger than 1 GeV say, the coupling 
constant at Q s will be small enough for small coupling methods to become applicable. In spite of this, even 
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. However, numerical solutions to the BK equation have been obtained in all 
[Ig , I20I I2H l22j . For a review of the properties of the numerical solutions cf. 



cases only [8j, [9|, |10|, 
momentum regions 

In the absence of exact solutions of the nonlinear evolution equations, the dipole cross section has been mod- 
eled using both theoretical considerations and experimental data. In the saturation region, such models have 
been derived theoretically on the basis of multiple rescattering of the dipole off the target, to various levels of 
sophistication. In what is now commonly referred to as the Color Glass Condensate (CGC) formalism [24], HH], 
the small- a; gluons are described as a color field generated by the large- a; partons, which act as static sources. 
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Scattering amplitudes are expressed as two-point functions of lightlike Wilson lines that are averaged over the 
possible source configurations, which evolve with x. The averaged two-point functions satisfy a tower of nonlin- 
ear evolution equations, the Balitsky-JIMWLK equations [26| , coupling them to ever higher n-point functions. 
However, in a mean field approximation the Balitsky-JIMWLK equations reduce to a single evolution equation, 
the BK equation. The saturated state of high-gluon density at small x described by the CGC has been under 
very active theoretical and experimental investigation in recent years. Valuable further tests will be provided by 
future colliders that are able to probe well into the saturation region, such as the LHC and a possible electron-ion 
collider. 

The phenomenological models for the dipole cross section that are used to investigate the available HERA and 
RHIC data are constructed on the basis of general theoretical arguments and contain only very few free parameters. 
In this way satisfactory descriptions of several types of processes have been obtained [27l . l28l . HiJ [3(| [3l|, [HJ . The 
data thus seem to confirm some general expectations concerning the dipole cross section, in and above the 
saturation region. Here we want to investigate whether the phenomenological models are compatible with the BK 
equation. We will compare the numerical solutions of the BK equation [22| with a recent model for the dipole 
cross section [29| that is intended to describe the momentum and x-dependence in a particular region above the 
saturation scale -the extended geometric scaling region. Our conclusions will also apply to the dipole models of 
[27l . I32I ] in that region. Regarding the cc-dependence it should be mentioned that the leading order BK equation 
is known [H, HH to lead to a faster evolution in x (Q 2 (x) ~ l/x x where A ~ 0.9 1 ) than the experimental data 
seem to favor (A ~ 0.3). This discrepancy can be reduced by introducing a running coupling constant and indeed 
in that case the DIS data seems to be describable taking both LO BK evolution and, for larger Q 2 , DGLAP 
evolution into account [35]. However, the consistency of including running coupling in an equation that is of 
leading order in a s is questionable. This indicates that going beyond the leading order BK equation may be 
required in order to arrive at a proper quantitative comparison. Here however, we will be mainly interested in 
the question of whether the qualitative features of the models are reflected by the solutions to the leading order 
BK equation, especially around the saturation scale and at subasymptotic rapidities. 

The outline of this paper is as follows. In Section II we discuss the properties of the phenomenological models 
of the dipole cross section in detail. In Section III the BK equation is briefly reviewed. Section IV is devoted to 
the study of the numerical solutions of the BK equation, by means of a comparison with a general Ansatz for the 
dipole cross section in terms of an anomalous dimension, in particular, the one proposed in Ref. !29|. We also 
study the ^-dependence of the saturation scale with and without running coupling and consider the dependence 
of the results on the initial conditions. In Section V we summarize the main conclusions that can be drawn from 
the comparisons we have performed. 



II. PHENOMENOLOGICAL STUDIES OF THE DIPOLE CROSS SECTION 



A very successful phenomenological study of experimental data using a model for the dipole cross sections was 
performed by Golec-Biernat and Wiisthoff (GBW) [3(|. The HERA data on the DIS structure function F 2 at 
low x (x < 0.01) they analyzed 2 could be described well by a dipole cross section of the form a — uoNcBwiTt:^), 
where the scattering amplitude Nqbw is given by 



V 



GBW 



(r t ,x) = 1 - exp 



-~/ t Ql{x) 



(i) 



r t denotes the transverse size of the dipole, cto — 23 mb and the x-dependence of the saturation scale is given by 

( A/2 

X 



.(.,-IG.V 



(2) 



with xq ~ 3 x 10 4 and A ~ 0.3. Although here we discuss DIS off the proton (A — 1), we mention that for 
nuclear targets Q 2 S increases as A 1 ^ . 



1 The exact value of A depends on a a = N c /na 3 , A ~ 0.9 corresponding to a 3 ~ 0.2. 

2 A further test of the dipole picture can be obtained from ratios of structure functions for which a bound has been derived recently 

M. 
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The above scattering amplitude displays a characteristic feature called geometric scaling. This means that it 
depends on x and r\ (r t being the Fourier conjugate of kt = Q) through the combination r\ Q 2 s {x) only. A model 
independent analysis of the DIS data from HERA shows that the low-.r data display geometric scaling for all 
Q 2 [13 , even though the GBW model, H])-([l]), was found to be inconsistent with newer, more accurate data [3lJ 
at large Q 2 . In Ref. [3l[ a modification of the GBW model was proposed which includes DGLAP evolution, as 
required to fit the Q 2 > 20 GeV 2 data. Alternatively, in Ref. [35[ the dipole scattering amplitude from the GBW 
model was replaced by a numerical solution to the leading order BK equation with running coupling and with the 
addition of a correction that satisfies DGLAP evolution in order to also describe the short distance behavior of 
N correctly. As emphasized in Ref. [32J (IIM), the solution to the BK equation (with running coupling) by itself 
only provides a satisfactory fit for relatively low Q 2 , up to a few GeV 2 . 

As mentioned, the DIS data show geometric scaling for x < 0.01 and all Q 2 (this was recently confirmed in 
Ref. [3q|). but the data above and below Q 2 /Q 2 (x) « 1 behave differently as a function of Q 2 /Q 2 {x) [13- On 
the basis of theoretical considerations one expects different geometric scaling behavior at low and high Q 2 /Q 2 
(a Gaussian-like function of Q 2 /Q 2 at low Q 2 /Q 2 and a power law fall-off at high Q 2 /Q 2 ). Moreover one 
expects an intermediate region, commonly referred to as the extended geometric scaling (EGS) region [3^.[39l[40|. 
starting from Q s , extending to a scale Q gs , where the geometric scaling behavior of the saturation region holds 
approximately. This EGS region grows with decreasing x. 

An estimate of Q gs has been obtained from the solution of the LO-BFKL evolution equation expanded to second 
order around the saturation saddle point po| . Here Q gs is taken to be the scale at which the violation of the 
geometric scaling behavior that holds at Q s becomes of the same order as the leading order scaling contribution. 
One then finds Q gs (x) ~ Ql(x)/A, where A is a nonperturbative scale of order Aqcd- A similar result for Q gs can 
be obtained from estimating the transition point between the LLA and DLA saddle points (Hj]. We emphasize 
that at Q gs one does not necessarily observe large deviations from geometric scaling, one may also reach a region 
with a different geometric scaling behavior. 

The theoretical studies of the EGS region have led IIM [HJ to propose a modification to the GBW model so 
as to obtain a satisfactory fit to DIS data in both the saturation region and the EGS region, without including 
DGLAP evolution at large Q 2 . In the same spirit, KKT [13 have considered a modification of the GBW model 
in the EGS region in order to obtain a good description of the pt spectra of hadron production at RHIC in d-Au 
collisions (for a review of saturation physics and d-Au collisions cf. [il|). Here we will restrict the discussion to 
the later modifications of the KKT model made by DH J [29( , but our conclusions for the EGS region will apply 
to the IIM and KKT models as well. We note that the saturation scale of the IIM model is to be multiplied by 
two in order to compare with the results obtained below. Other dipole models are discussed in Refs. [42l. l43l|. 

The model proposed by DHJ offers a good description of the pt spectra of hadron production at RHIC in d-Au 
collisions in both the midrapidity and forward regions, and even, as has been shown recently [44j . the p-p data in 
the very forward rapidity region. 

The modified dipole scattering amplitude in the EGS region is given by [13, HI, : 



N(r t , x) = 1 — exp 



l(r 2 Q 2 (x)r irt - x) 



(3) 



The exponent 7 is usually referred to as the "anomalous dimension" , although the connection of N with the gluon 
distribution may not be clear for all cases considered below. This anomalous dimension 7 changes the evolution 
in x, thus governing the possible violation of geometric scaling for scales between Q s and Q gs . The unitarized 
form of N is important only for scales around and below Q s . For larger scales one can use to good approximation 

N{r t ,x)^\{r 2 t Q 2 s {x))^\ (4) 

which allows one to make a connection with the DGLAP region 3 , where iVDGLAp(fct, %) °^ (N being the 
Fourier transform defined in Eq. ([?]))■ which is the tail behavior determined by one gluon (i-channel) exchange. 
Hence, the "DGLAP" limit of 7 would be 7 -> 1. 

DHJ (following in part IIM and KKT) impose the following requirements to determine a parameterization of 
7 in the EGS region: 



3 More specifically, at low x and high Q 2 a logarithmic evolution in both these variables may be required. 
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1. For any fixed x, in the limit r t — > 7 should approach 1. 

2. At the saturation scale, 7 should be constant, so that N displays exact geometric scaling: ^(r t = l/Q s , x) 
7s. The constant 7 S is chosen to be ~ 0.628. 



3. If one writes 7 = 7 S + A7, then A7 should decrease as 1 jy for y — > 00 at fixed r 2 Q 2 . 

4. At large, but fixed, rapidity y, the extent of the geometric scaling window should be consistent with the 
estimated scale Q gs = Qs/Aqcd- 

A few comments on these requirements are in order. Requirement 1 is imposed in order to reproduce the "DGLAP" 
limit N ~ r\ . Requirement 2 and 3 were motivated by properties of the BFKL equation [3^]. The idea is that 
since in the saddle point approximation the solution of the BFKL equation exhibits the property of geometric 
scaling, one can use it as an approximation to the solution of the BK equation near the saturation scale Q Sl 
provided that one imposes an appropriate boundary condition at Q s |33l . [34l . |40| . commonly referred to as a 
saturation boundary condition. Requirement 4 [23 ] is of a more phenomenological nature, and determines the 
speed with which 7 increases from 7 S ~ 0.628 at Q s to a value near 1 at Q gs w Qs/Aqcd (which is a function of 
y through Q s )- As can be seen from e.g. Ref. |4lj |. requirement 3 is a prerequisite for requirement 4. 
Following the above requirements, the anomalous dimension 7(r t , x) of DHJ is parameterized as (29| 

< \ \og{l/r 2 Q 2 s (x)) 
j{r t ,x) = 7s + (l-7 s ) — - 5 
Ay + d^y + \og(l/r z t Qj(x)) 

where y = \ogx$/x is minus the rapidity of the target parton. The saturation scale Q s {x) and the parameter A 
are taken from the GBW model, as given in Eq. Requirement 4 is satisfied, which can be seen as follows: 
for large, but fixed rapidity (y ^> (d/X) 2 ), the EGS region roughly extends to \ogQ 2 s /Q 2 ~ Ay, which implies 
that numerically Q gs w Q 2 /Aqqd- The parameterization ([5| includes subleading corrections ~ d^fy which govern 
geometric scaling violations at subasymptotic rapidities. For d ~ 1.2 a good descri ptio n of d-Au data from RHIC 
over a considerable rapidity range, from mid- to forward rapidities, was obtained [291 ] - This parameter will play 
no role in our discussion, as we do not aim to obtain the best possible expression for either N or 7, but rather we 
want to investigate whether the above type of parameterization is compatible with the BK equation. 

Two further aspects of DHJ's Ansatz ^ for N with the anomalous dimension of Eq. should be mentioned. 
First, in order to calculate observables (and to connect to the BK equation), one needs to formulate the dipole 
scattering amplitude in momentum space by performing a Fourier transform. To simplify this procedure, DHJ 
(following a similar step by KKT) replaced 7(7-4,2;) by ~f(l/kt,x). This replacement should be viewed as an 
approximation that becomes better as k increases or equivalently, as r decreases. This simplifying approxima- 
tion could be regarded as part of the Ansatz. Second, quarks and gluons are described by different scattering 
amplitudes, since they couple differently to the gluons. The above Ansatz ^ was applied by DHJ to gluons 
and is referred to as Na- The corresponding expression for quarks is called Np and is obtained from Na by the 
replacement Q 2 — > (Cf/Ca)Q 2 s = (4/9)Q^. In the comparison to the GBW model, which is a model for Np, one 
should take into account this factor, leading to Q s (xo) = 1 GeV instead of 0.67 GeV. 

Eq. was intended to parameterize 7 in the geometric scaling window only, so the question arises what forms 
N takes in the saturation region, where r t > l/Q s - For instance, it has been suggested [44( to be of the form ([3]) 
with 7 a constant, equal to y s , so that the Ansatz for N displays exact geometric scaling throughout the saturation 
region. However, this idea has the drawback that the resulting N is not a smooth function of r t , a problem also 
present in the IIM model. Its Fourier transform would therefore show artificial oscillations. To resolve this issue, it 
seems better to use the BK equation (or other, more appropriate, evolution equations) to determine the behavior 
of N in the saturation region. The goal of this paper is to study in more detail the relation between the DHJ 
Ansatz and the BK equation. We use the "BKsolver" program [45( to obtain a numerical solution to the BK 
equation and study whether the dipole profile of DHJ and similar models is consistent with this solution. Also, 
we will address the behavior of N in the saturation region as dictated by the BK equation. More specifically, 
we will investigate whether the requirements 2 and 3 on 7, which were motivated by the BFKL equation with a 
saturation boundary condition, are consistent with the nonlinear evolution of the BK equation. 
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III. BALITSKY-KOVCHEGOV EQUATION 

The BK equation for N (more specifically, for Np) reads @, 0] 



dN((x t ~y t ) 2 ,x) a s f 2 (x t - y t f 



d'z t - 



N{{x t -z t ) 2 ,x)+N{(z t ~y t ) 2 ,x) 



dy 27r J (x t - z t ) 2 (y t - z t ) 2 

- TV {{x t - y t ) 2 ^) " ^ {{S t - z t ) 2 ,x) N {{^ - y t f,x) ] . (6) 

Here a s = a s N c /i:. We will not consider the impact parameter dependence of N in this paper. 
Two types of Fourier transform are often considered in this context: 

N{k,x) = f ^e rk ^N{r,x), (7) 



Af{k,x) = [ ^e**^, (8) 



rf and k 2 — fc 2 . The first (regular) Fourier transform is the one that enters the cross section 
descriptions, such as the one considered by DHJ [28|, [2i|. The second transform, which includes the additional 
factor 1/r 2 , is more directly related to the (unintegrated) gluon distribution, as explained in detail in Ref. [5l| . 
In terms of Af, the BK equation can be very compactly written as 

d Y M = x {-d L )M - N 2 , (9) 

where Y = ya s , L = log(fc 2 /fco), f° r some arbitrary scale fco, and xil) IS the BFKL kernel, 

X(7) = 2^(1) -#y)~ ^(1-7)- (10) 

In the limit of large k Q s ), or equivalently, small r, the nonlinear BK equation © reduces to the linear 
BFKL equation. In the saddle point approximation, the solution of the BFKL equation behaves as N{r, x) ~ r 27 
for small r, hence the exponent 7' determines the small-r limit of the anomalous dimension 7 through Eq. f¥|). 
Therefore, 7 is determined by either the BFKL saddle point 'Js.p. or the initial condition: if in the limit r — > the 
starting distribution has the property N{r, Xq) ~ r 270 , then for any fixed x the asymptotically small r behavior 
either remains r 2la if 70 < Js.p., or becomes r 27a p - if 70 > 7s. p.- We will demonstrate this observation explicitly 
below. For finite x, the BFKL saddle point approaches j s . p . w 1 in the limit of small r, cf. e.g. [40j. Therefore, 
any solution to the BK equation will only satisfy requirement 1 of DHJ if the starting distribution in the small-r 
limit has the property that 70 > 1. 

Around the saturation scale, the solution of the BK equation has been approximated by a solution of the 
BFKL equation which has been required to satisfy a saturation boundary condition [H, H3, [4C| ■ This leads to the 
aforementioned saddle point at saturation, 7 S ps 0.628, which is the solution of the equation x{ls)/ls — x'ils) 
[i],|4(|. However, unlike in the small-r limit where the BFKL solution and Eq. §S§ reduce to the same form, around 
the saturation scale there is no clear correspondence between this saddle point and the anomalous dimension of 
Eq. ©. 

Alternatively, one can first approximate the BK equation ^ by expanding the kernel x around the saddle 
point, and subsequently solve the resulting equation analytically. This has been done in Refs. [§, [l(| HH, [3 where 
a traveling wave solution has been obtained that is valid in a finite range of fc, which grows with the rapidity. 
The value 7 S w 0.628 enters the analytic expressions for Af{k,x) of the traveling wave solutions 4 . However, also 
for this solution 7 S does not correspond directly to the anomalous dimension of Eq. at the saturation scale, 
since there is no direct relation between N{r, x) and M{k, x) at the respective points r = 1/Q S and fc = Q s . 

We conclude that requirement 2 of DHJ, which sets 7 at the saturation scale equal to 7 S » 0.628, does not 
follow from these analytical results. For completeness we mention that other analytical aspects of the solutions 
of the BK equation, such as the small-fc limit, have been investigated as well pi, Ha, [T3| . 



4 For initial conditions with 70 < 7s , 7s is replaced with 70 in the traveling wave expression. 
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Here, we will restrict ourselves to comparing the phenomenological DHJ model with numerical solutions. Such 
solutions have been obtained by several groups [H, Ufl E3, Gil EH 0, l2ll l22l ]- even including impact parameter 
dependence [13, IH, [4!| [5(| . The code employed by the authors of Ref . [22[ to numerically solve the BK equation 
has been made publicly available [45| . We will use this code in our subsequent investigation of the DHJ Ansatz 
and its possible extension to the saturation region. But first we want to point out that these numerical solutions 
do not show geometric scaling for finite rapidities. The deviations can become significant even in the saturation 
region at moderate y. As is known, in the limit of large rapidities, geometric scaling is recovered; M(k, x) is 
given by a function Moo that depends on k/Q s (x) only. Fig. [T^ shows the relative difference between M and the 
asymptotic solution Moo as a function of k/Q s {x). The saturation scale Q s (x) is defined as the value of k where 
M equals some arbitrary fixed value for which M 2 becomes non-negligible. Based on our study in Sec. IV. B we 
have chosen M(Q s (x),x) = 0.2. Due to this definition M — Moo vanishes at the saturation scale. Clearly, the 
scaling violations in the saturation region are considerable for finite rapidities. Fig.[T|D shows the dipole scattering 
amplitude N(r,x) in coordinate space as it follows by a Fourier transformation (fTTI) from M(k,x). As expected, 
N(r, x) vanishes for small r and approaches 1 in the saturation region r 3> 1/Q s (x). Like in momentum space 
N(r, x) approaches a limiting curve that depends only on rQ s (x) and geometric scaling is recovered. Throughout 
this paper we have chosen a constant coupling a s = N c /Tra s ~ 0.2. 




FIG. 1: a) The relative difference 1 — Moo(k/Q s {x))/M(k/Q s (x),x) between the function M(k,x) © following from 
the BK equation and its geometrically scaling asymptotic lim^^o M(k/Q s (x), x) — Moo{k/Q 3 (x)) for various rapidities 
y = \ogx /x. 

b) The dipole scattering amplitude N(r, x) resulting from a Fourier transformation (|lip of M(k, x) as a function of rQ s (x) 
for various rapidities y — \ogxo/x. 



IV. ANOMALOUS DIMENSION FROM THE BALITSKY-KOVCHEGOV EQUATION 
A. Anomalous dimension as a function of r and x 



The BKsolver program [4o| provides numerical values for the amplitude M(k, x) as a function in momentum 
space. In order to use this solution of the BK equation to constrain j(r, x), one first has to find N(r, x) by Fourier 
transforming back to coordinate space: 



N(r,x) = r 1 



d 2 k t 
2tt 



M(k,x)=r 2 dkk J (kr)M(k,x) . 



Using the Ansatz §S§ one can extract j(r, x) from the resulting N(r, 



7(r, x) 



log 



log 



(1-N(r,x)) 4 



log[r2Q2(x)] 



(11) 



(12) 
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This equation requires as a separate input the value of Q s (x), which can be found as follows. One can fix Q s (x) 
in such a way that at the saturation scale the amplitude obtained from the BK equation (fT2"|) equals the Ansatz 
©, which at Q s becomes independent of 7, 

N(r = l/Q a , x)=l- exp [-1/4] w 0.22 . (13) 

Equating this fixed value with the general relation (fTTjl allows one to calculate the saturation scale Q s (x) by 
solving the following equation, 

1 f°° 

— — / dkkJ (k/Q s (x))M(k,x) = l-exp[-l/4}. (14) 

Qsix) Jo 

Combining the resulting values of the saturation scale with Eq. (|12[) . we obtain a numerical result for 7(7*, x), 
which is shown in Fig. [2] as a function of l/(rQ s (x)). This means that the saturation region is located to the left 
of 1 on the horizontal axis, so that it is easier to compare with later figures that display related quantities as a 
function of k/Q s . Note that although the Ansatz (J3j) was specifically intended to describe N in the EGS region, 
we also show 7 in the saturation region. 



y=o 

y=l 




0.3 1 1 ' ' 

0.1 1 10 100 1000 

V(rQ.(ar)) 

FIG. 2: >y(r,x) resulting from the relations JTT} and (|12|l as a function of l/(rQ s (x)) and y = logzo/a:. 

The resulting 7(7-, x) has the following features: 

1. For r — > 0, 7(7", x) asymptotically approaches 1. 

2. At the saturation scale, 7(7*, x) is not a constant. 

3. For decreasing x, 7(7", x) approaches a limiting curve, indicated in Fig. [2] by y = 00. Hence, after a longer 
evolution one indeed recovers geometric scaling. 

The fact that for small distances 7 asymptotically approaches 1 is understandable from the BK equation, since 
in this limit it reduces to the BFKL equation. As mentioned before, in the limit of small distances, the solution 
to the BFKL equation is dominated by either the saddle point or the initial condition, both leading to 7 — * 1, 
since here we use the MV model ([17]) as the initial condition. 

Requirement 2 of the DHJ Ansatz is clearly not satisfied by the numerical solution of the BK equation at 
the saturation scale. Instead of writing 7(7*, x) — "/ s + &y(r,x), it seems more natural to consider the following 
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split-up: 



j(r,x) = ~foo(rQ s (x)) + Aj(r,x) , 

where it turns out that, similar to requirement 3, A7(r, x) decreases as 1/y for y 
saturation scale 7 is given in the small- a: limit by, 

lim 7 (r - 1/Q s {x), x) = 7oo (l) « 0.44 , 

x— >0 



(15) 

00 and fixed rQ s {x). At the 



(16) 



a value that is significantly below 7s — 0.628. 

It would be interesting to see whether the numerically obtained 7 could fit the available data as well, which 
does not seem unlikely since KKT employed 7s = 0.5 and were able to fit d-Au data, but this is not the aim of 
the present paper, where we instead focus on the comparison of existing models for the dipole cross section with 
the numerical solution of the BK equation. Moreover, a word of caution has to be added about the results for 
small values of y (below y = 5, say): the numerical solutions are obtained after a choice of a starting distribution 
at y — (x — xq). For "short" evolutions (small y- values), the properties of the starting distribution are still 
visible in the result, i.e. not only in the small-r limit. The starting distribution should be appropriate to the 
choice of xo, which itself should be small enough for BK evolution to be applicable. If xo is sufficiently small, 
or alternatively if A is sufficiently large, one can take for instance the McLerran-Venugopalan (MV) model [24| . 
One could also consider taking the DHJ model, which however is restricted to the EGS region, therefore requiring 
some extrapolation to the entire r or k range. For larger values of xq one may consider using one of the standard 
parton distributions, although these are not defined over the whole range either. 

Due to this problem, throughout this paper we have used the MV model as a starting distribution 5 : 



N(r, x) = 1 — exp 



1 



(r^(aO)tog(e+l/(r'A J )) 



(17) 



which is one of the distributions considered in Ref. [22j. In Sec IV. E we show some results for another choice of 
starting distribution, to illustrate the influence of the initial condition. The dependence of the numerical results 
on the starting distribution has also been investigated in Ref. 22], to which we refer for additional information. 

Hence, for small rapidities (y < 5) one should avoid drawing too strong conclusions from a comparison of the DHJ 
parameterization with the numerical results that were obtained with the MV starting distribution. Nevertheless, 
we will be able to draw some qualitative conclusions also for lower y- values from our results in Sec IV. B. 



B. The replacement y(r, x) r y(l/k,x) 

As mentioned before, DHJ actually considered Eq. §S§ with j(r,x) replaced by 7 (l/fc, x) [13, H^, [29j. This 
approximation scheme we will discuss next. The procedure of extracting 7 becomes quite different when 7 
depends on k, since the dipole cross section N then depends on both r and k, so that it is not related to Af(k, x) 
by a straightforward inverse Fourier transform pip anymore: 



Af{k,x) 



dr 

— J {kr) 
r 



1 — exp 



-Vtf(*)) 7(M) 



(18) 



Instead of by the inverse Fourier transform, we will extract 7 by trying to numerically solve Eq. (|18p. imposing 
the following condition. In order to test the Ansatz of DHJ, we will fix ^{k, x) in such a way that it equals the 
constant 7 S at the saturation scale: 



7(Q»(a:),ar) = 7s ~ 0.628. 



(19) 



5 There are many ways in which the infrared divergence that occurs in the MV model can be regularized; the replacement of 
log(l/(r 2 A 2 )) — * log(e + l/(r 2 A 2 )) is one option. In Ref. it is shown that the results do not depend much on the choice of 
regularization. 
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This determines the dependence of Q s on x, which can be found by explicitly solving 



Af(Q s (x),x) = I — JqO) ( 1 - exp 



4 V ' 



0.19. 



(20) 



Now we can extract 7 from relation (fl8|) for any given value of x and k. Fig [3] shows the results for j(k, x) as a 
function of k/Q s in the EGS region, for a broad range of rapidities. For small rapidities the resulting 7 looks very 
similar to the one of DHJ (cf. Fig. 4 of Ref. [29J). This gives us an indication of how the DHJ distribution would 
behave under BK evolution. We reiterate that one cannot use the DHJ distribution itself as a starting distribution 
as it is not defined over the whole r or k range. As one can see, for larger y the resulting 7 is not compatible with 
the DHJ parameterization ([5]) anymore; it first decreases (A7 < 0) before it rises towards 1 asymptotically. 




0.5 1 1 1 1 1 1 1 1 

15 10 15 20 

k/Q s (x) 

FIG. 3: j(k,x) as a function of k/Q s (x) for various rapidities y = \ogxo/x. 

In Fig. H] the same result is shown but now for a smaller choice of constant 7 at the saturation scale, 
7(fc = Q s ,x) — 0.44. This value is the small- a; limit (fT6|) of j(r = 1/Q s (x),x) in the approach of the previ- 
ous subsection. In this case A7 remains positive above the saturation scale and the result looks very much like 
the DHJ parameterization for all y. A fit of 7 as given in Eq. (f5]) to the numerical results for k/Q s = 1, 2, . . . , 5 
and rapidities y = 1, 2, . . . , 5 (in order to allow for a comparison to DHJ's results) yields Awl and d w 3, although 
it must be emphasized that the shapes of the curves are not exactly of the form (|5|). 



C. The saturation region 



The DHJ Ansatz was specifically intended to describe the dipole scattering amplitude in the EGS region, but it 
is nevertheless interesting to see what 7 will look like if one were to assume the same form of N in the saturation 
region. One may expect the solution to be geometrically scaling in the region r > 1/Q S , which would be the case 
if 7(7", x) — r y(rQ s (x)). A constant 7 would of course be a special case. However, the numerical solution of the 
BK equation does not display geometric scaling, not even in the saturation region, although for large rapidities 
Af(k,x), as a function of k/Q s (x), converges to a well-defined limit. The same conclusion applies to 7 upon 
assuming the Ansatz ^ for N(r,x). This is clearly seen in Fig. [3J 
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k/Q s (x) 

FIG. 4: y(k,x) as a function of k/Q a (x) for various rapidities y = \ogxo/x. 

The situation is similar if we replace j(r,x) with j(l/k,x), keeping 7 at the saturation scale k — Q s fixed, with 
the exception that in some cases no solutions for 7 can be found at all in the saturation region. The existence of 
a solution in a certain range of k/Q s depends on the chosen saturation scale or, in the above discussed approach, 
equivalently, on the fixed value of 7 at the saturation scale, which determines the saturation scale. In fact, we 
find that if j(k = Q s (x), x) is fixed to equal 7 S = 0.628, the Ansatz (J3j) is not compatible with the solution of the 
BK equation in the range k/Q s « 0.2 — 0.4, where no value of ^(k, x) satisfies the relation (fl"8"|) for large rapidities. 
If instead we fix 7 at the saturation scale to equal the limit (fTI))) , the saturation scale is shifted towards smaller 
values, so that there is a unique solution for 7 for all k and x. Above the saturation scale 7(fc,x) is in agreement 
with the type of parameterization of DHJ (|5|). For asymptotic rapidities, 7 approaches a limiting curve, which is 
geometrically scaling for all k values, not just at the saturation scale (which was the case by construction). Note 
that below the saturation scale the solution is far from constant, although for large rapidities the solution is rather 
constant for a considerable momentum range below Q s . But eventually, for very small k, it drops significantly to 
a small nonvanishing value. We emphasize that this is the case under the assumption that the form §3§ also holds 
for momenta below the saturation scale. Let us mention that in the region where 7 eventually drops towards 
smaller values r Q s (x) becomes effectively already so large that N(r, w 1, so that the exact value of 7 eventually 
becomes irrelevant, see Fig. [T|d. 

D. Saturation scale dependence on x 

In Fig. [5]we display the values of Q s {x) for a fixed coupling a s — 0.2 as a function of y = logxo/x, in order to 
show that (at least for larger rapidities, y > 10 say) the behavior Q s (x) ~ x~ x / 2 is found in all scenarios considered 
here. The values of the parameters can be compared with the GBW values of A ~ 0.3 and Q s (y — 0) ~ 1 GeV. 
We note that strictly speaking one has to include a factor 2/3 in Q s (y = 0) if the DHJ Ansatz for Np is used for 
calculating Q s {x). This is the case in e.g. Eqs. (fl4")) and (|2H)) . The larger value of A » 0.9 is in agreement with 
the asymptotic result A ~ a s x{^ s )/ls ~ a s 4.88 from LO BFKL equation with saturation boundary conditions 
[11] and A = 4a s from the DLA approach [ll, HH, [H, . 

We have also investigated the case where the coupling constant runs with the energy scale, which is given by the 
inverse of the relevant dipole size. For details, we refer to [I2,[45j]. The energy scale turns out to be effectively of 
the order of the resulting saturation scale Q s . The running of the coupling results in a change in the dependence 
of Q s (x) on x, for which we now find approximately Q s (y) = Qs(0) exp[c»/y/2], which in fact is equally well 
supported by all relevant DIS data [38]. Since we have plotted every result for the anomalous dimensions as 
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- A = 0.90,Q S (0) 


= 0.34 GeV 




A = 0.91,0,(0) 


= 0.42 GoV 




A = 0.91,(3,(0) 


= 0.29 GoV 
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QAv), Fig- 2 




□ 


Q,(y), Fig. 3 
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Q4v), Fig. 4 





30 

V = logz Q /x 



50 



FIG. 5: The saturation scale from the various approaches for y = \ogxo/x = 0, 5, 10, 
of the form Q s {y) = Qs(0) exp[A y/2] in the fixed coupling case, where a 3 — 0.2. 



. , 50 (symbols) and the related fits 



a function of dimcnsionless quantities, i.e. rQ s (x) or k/Q s (x), none of those results is significantly affected by 
including running in the coupling constant. It has to be mentioned that since the BK equation is of leading order 
in a s , including the running coupling is of course not a consistent treatment of higher order effects. For that 
reason we mainly focus on the fixed coupling case. 



E. Dependence on the initial conditions 



To illustrate the dependence of some of our results on the choice of a starting distribution, we replace the MV 
distribution with the following choice: 



N(r, x = xq) = 1 — exp 



(21) 



Fig. shows the resulting 7(7*, x) as a function of \/(rQ s (x)) for the initial conditions 70 = 0.6,1 and 1.1, 
in a very large r range in order to display the asymptotic small-r behavior. As expected from the discussion in 
Sec. HID in the limit of very large l/(rQ s ) at finite rapidity, 7 approaches either the initial condition 70, if 70 < 1, 
or the saddle point, which is 1 in this limit. 

Fig. [lb shows a similar plot for the initial conditions 70 = 0.5,0.6 and 0.7 at larger rapidities. For 70 = 0.5 in 
the large- y limit 7 approaches a limiting function that monotonically rises towards 70 at small rQ s . The limiting 
function for 70 = 0.6 behaves similarly. For 70 — 0.7 the limiting function coincides to good approximation with 
the one of 70 = 0.6 up to small values of rQ s , but becomes slightly larger as rQ s decreases further. In the 
entire investigated range the limiting function of 70 = 0.7 coincides with the one of the MV starting distribution 
(|17p. partly shown in Fig. [5J We conclude that if ^o( r Qs) is sufficiently large, one always reaches a universal 
limiting function "fcoi^Qs)- This function equals approximately 0.44 at the saturation scale and approaches 
a limit 7oo(rQ s — > 0) which seems to be slightly larger than 0.6. As discussed below, the latter limit seems 
consistent with the theoretical expectation of 0.628. In the specific case of constant 70 the universal limiting 
curve is always reached if 70 is larger than ^ooi^Qs — > ). 

A similar dependence on 70 has been observed in Ref. [l0( on the traveling wave solution. This solution depends 
on a parameter 7 which equals either 7^ or 70, whichever is the smallest. In our approach the role of "f s of the 
traveling wave is played by the function 7^. At large k/Q s and y the traveling wave solution is dominated by 
the term (k 2 /Q^(x))" 1 , leading to the previously mentioned expectation that 700 (fQs — ► 0) = 7 S . 
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FIG. 6: a) 7(r, x) for three different rapidities and initial conditions 70 = 0.6, 1, 1.1. To illustrate the small r asymptotics, 
we show 7(7-, x) over a very large range. Note that the curves with 70 = 1 and 1.1 for y = 50 are hard to distinguish, 
b) j(r, x) as a function of l/(rQ s (a;)) and initial conditions 70 = 0.5,0.6 and 0.7 for large rapidities y = logxo/x. Here, 
the curves for y — 60 and y = 105 largely coincide for both 70 = 0.5 and 0.6. 



Finally, we compare our results to those of Ref. [4§], especially Fig. 4. There it is found that for small rQ s 
the scattering amplitude can be described by N(r) = a(rQ s ) 27 (log(rQ s ) + (5), where 7 approaches approximately 
0.65 for large y, regardless of the initial condition. We find that the particular value of 0.65 is consistent with 
our result for 700, if we take into account the different functional form for N(r, x) and the specific range of r and 
y in which 7 was calculated. However, even for 70 = c/2 = 0.42, Ref. 14911 find no dependence of 7 on the initial 
condition, which seems to disagree with our results and those of Ref. • 



V. CONCLUSIONS 



The numerical solutions of the BK equation do not display exact geometric scaling, although they approach 
a solution showing such scaling at asymptotic y. Assuming the solutions to be of the form Q, where scaling 
violations are encoded in the "anomalous dimension" 7, therefore leads to the conclusion that 7(7*, x) is not a 
function of rQ s (x) exclusively. In particular, it is never simply a constant, not even at the saturation scale 
(r = l/Q s ). At asymptotically large rapidities, 7 reaches a limiting function 7oo(rQs(x))- This function is 
universal for a large range of initial conditions. At the saturation scale, this function equals approximately 0.44, 
which is considerably smaller than the corresponding values of 7 in the phenomenological models [13, HH, H3] . 
For small values of rQ s the limiting function seems to reach j s , in accordance with the traveling wave results of 
Refs. 

This is not in contradiction with the behavior of 7 in the limit of rQ s (x) — > at any fixed x, where 
7 is determined by the starting distribution or the BFKL saddle point. It merely implies that the two limits 
rQ s (x) — > and x — > do not commute. 

Although the numerical solutions of the BK equation lead to a 7(7", x) that is never a constant as a function of x at 
the saturation scale, performing the replacement of 7(7*, x) — > 7(l/fe, x) does allow one to find a solution for which 
7(fe = Q s ,x) is kept fixed. The behavior of 7(1/^,2:) is then qualitatively similar to the DHJ parameterization 
for small rapidities. However, the usually considered choice "f(k = Q s ,x) = "f s = 0.628 yields some unwanted 
features, i.e. the fact that A7 < in a region above the saturation scale and the absence of solutions below the 
saturation scale, although the Ansatz was not intended for that region. Keeping j(k = Q s ,x) fixed at a smaller 
value, e.g. at joc(rQ s = 1) w 0.44, seems more suitable, but it remains to be investigated whether such a choice 
allows for a good fit of all relevant DIS, d-Au and p-p data. One might expect this to be possible, as the KKT 
model [27|, which has 7 = 0.5 at the saturation scale, is able to describe the d-Au data. 
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The resulting saturation scales in the various approaches we adopted, evolve quite similarly with decreasing x, 
namely approximately as Q s (x) = Q s (xq){xq / x) x ^ 2 , albeit with somewhat different normalizations Q s {xq). In 
the case of a running coupling constant, the de pen dence on y changes from exp(Ay/2) to exp(c v /y/2), which is 
equally well supported by all relevant DIS data [381 ] - 




Our conclusions about the DHJ model apply as well to the IIM and KKT models in the EGS region. It would 
be interesting to consider modifications of these models for the dipole scattering amplitude that are compatible 
with both the BK equation and the data. Given the fact that the BK evolution does not respect geometric 
scaling around Q s , phenomenological parameterizations that reflect this feature would seem a natural choice. 
Fortunately, the LHC and a possible future electron-ion collider will provide data over a larger range of momenta 
and rapidities, so that one can expect to test the evolution properties of the models more accurately. 
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